Production of n^pK^ in p+p reactions at 3.5 GeV 

beam energy 

G. Agakishievf, A. Balanda^ D. Belver^ A. Belyaev^, J.C. Berger-Chen^ 
A. Blanco^, M. BohmerS J. L. Boyard°, P. Cabanelas'i, E. Castro"!, 
S. Chernenko^, T. Christ', M. Destefanisj , F. Dohrmann'^, A. Dybczak"^, 
E. Epple"^, L. Fabbietti'^, O. Fateev^, P. Finocchiaro'^, P. Fonte'^''^, J. Friese', 

I. Frohlich^, T. Galatyuk^'*, J. A. Garzon^, R. Gernhauser', C. Gilardij, 
M. Golubeva\ D. Gonzalez-Di'az", F. Guber\ M. Gumberidze°, T. Heinz'^, 
T. Hennino°, R. Holzmann'^, A. lerusalimov^, 1. lori'^'^, A. Ivashkin', 
M. Jurkovic', B. Kampfer*^'^, K. Kanaki'^, T. Karavicheva\ I. Koenig*^, 
W. Koenig'^, B. W. Kolb^, R. Kotte^ A. KrasaP, F. KrizekP, R. Kriicken', 
H. Kuc"'°, W. Kiihnj, A. Kugler?, A. Kurepm\ R. Lalik^^, S. Lang*^, 
J. S. Langej, K. Lapidus^, T. Liu°, L. Lopes^, M. Lorenz^, L. Maier', 
A. Mangiarotti'^, J. Markert^, V. MetagJ, B. Michalska^ J. Michel^, 
E. Moriniere°, J. Mousa"^, C. Miintz^, R. Miinzer'^, L. Naumann'^, 
J. Otwinowski'^ , Y. C. Pachmayer^, M. Palka"^, Y. Parpottas'''™, 
V. Pechenov^, O. Pechenova^, J. Pietraszko^, W. Przygoda"^, B. Ramstein°, 
A. Reshetin\ A. Rustamov^, A. Sadovsky\ P. Salabura'^, A. Schmah'', 
E. Schwab'^, J. Siebenson'^, Yu.G. SobolevP, S. Spataro"", B. Spruck-i, 
H. Strdbeles, J. StrothS'<^, C. Sturm^, A. TarantolaS, K. Teilab§, P. TlustyP, 
M. Traxler<^, R. Trebacz^ H. Tsertos'^, V. WagnerP, M. Weber', 
C. Wendisch*^, J. Wiistenfeld*^, S. Yurevich*^, Y. Zanevsky^ 

""Istituto Nazionale di Fisica Nucleare - Laboratori Nazionali del Sud, 95125 Catania, 

Italy 

^ LIP-Laboratorio de Instrumentagdo e Fisica Experimental de Particulas , 
3004-516 Coimbra, Portugal 
Smoluchowski Institute of Physics, Jagiellonian University of Cracow, 30-059 Krakow, 

Poland 

CSI Helmholtzzentrum fiir Schwerionenforschung GmbH, 64291 Darmstadt, Cermany 
'^Institut fur Strahlenphysik, Helmholtz-Zentrum Dresden-Rossendorf, 01314 Dresden, 

Germany 

^ Joint Institute of Nuclear Research, 141980 Dubna, Russia 
^ Institut fiir Kernphysik, Goethe- Universitdt, 60438 Frankfurt, Germany 
^Excellence Cluster 'Origin and Structure of the Universe' , 85748 Garching, Germany 



*laura. fabbietti@ph.tuiii.de 
* * j ohannes . siebenson@ph. turn, de 



Preprint submitted to Elsevier 



February 14, 2012 



^Physik Department El 2, Technische Universitdt Miinchen, 85748 Garching, Germany 
Ul.Physikalisches Institut, Justus Liebig Universitdt Giessen, 35392 Giessen, Germany 

'^Istituto Nazionale di Fisica Nucleare, Sezione di Milano, 20133 Milano, Italy 
''Institute for Nuclear Research, Russian Academy of Science, 117312 Moscow, Russia 
"^Frederick University, 1036 Nicosia, Cyprus 
'^Department of Physics, University of Gyprus, 1678 Nicosia, Cyprus 
"Institut de Physique Nucleaire (UMR 8608), CNRS/IN2P3 - Universite Pans Sud, 

F-91406 Orsay Cedex, France 
P Nuclear Physics Institute, Academy of Sciences of Czech Republic, 25068 Rez, Czech 

Republic 

'^LabCAF. Dpto. Fisica de Particulas, Univ. de Santiago de Compostela, 15706 Santiago 

de Compostela, Spain 
^Also at Lawrence Berkeley National Laboratory, Berkeley, USA 
^Also at ISEG Coimbra, Goimbra, Portugal 
*Also at ExtreMe Matter Institute EMMI, 64291 Darmstadt, Germany 
'^Also at Technische Universitdt Darmstadt, Darmstadt, Germany 
'"Also at Technische Universitdt Dresden, 01062 Dresden, Germany 
'"'Also at Dipartimento di Fisica, Universitd di Milano, 20133 Milano, Italy 
^Also at Dipartimento di Fisica Generale and INFN, Universitd di Torino, 

10125 Torino, Italy 



Abstract 

We study the production of Ti'^n^pK'^ particle quartets in p+p reactions 
at 3.5 GeV kinetic beam energy. The data were taken with the HADES 
experiment at GSI. This report evaluates the contribution of resonances like 
A(1405), E(1385)°, A(1520), A(1232), N* and K*^ to the J:VpK+ final 
state. The resulting simulation model is compared to the experimental data 
in several angular distributions and it shows itself as suitable to evaluate the 
acceptance corrections properly. 
Keywords: A(1405), p+p collisions 
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1. Introduction 



For already half a century the A(1405) is a well known resonance with 
strangeness S = —1, Isospin 1 = and spin |. Even though its four star 
character suggests a good understanding of this baryon, its inner structure is 
still a topic of investigation. Indeed it is difficult to describe the A (1405) as 
a three quark baryon, as it is lighter than its nucleon partner, the A^*(1535). 
Also the large mass difference to the A(1520) can not be understood in terms 

n 

of spin-orbital coupling With the mass of the A (1405) lying slightly below 
the KN threshold another picture of this particle was established. From the 
analysis of the KN scattering length Dalitz and Tuan predicted the A(1405) 
in 1959 , already two years before its experimental discovery. Nowadays 
the A(1405) is described in a coupled channel approach based on chiral dy- 
namics J]. Here this baryon is generated dynamically as an interference of 
two states, a K~p bound state and a Svr resonance. However, this two pole 
structure cannot be observed directly in the Svr invariant mass spectrum, as 
the Svr pole is located far in the imaginary part of the complex energy plane. 
With these predictions, the structure of the A(1405) is interesting in terms 
of a deeper understanding o 
data are available for vr^+p 



)f the kaon-nucleon interaction. Experimental 

inn 

si], K +p [6| and 7+p [2] reactions. First re- 
sults on p+p data were reported in [8]. But only the neutral decay channel 



(A(1405) — )■ S°7r°) was investigated. We also study p+p reactions and con- 
centrate on the charged decay channels (A(1405) — )■ S^'^vr^). In order to 
extract precisely the spectral function of the A(1405), all the reactions that 
contribute to the T,'^n^pK~^ particle quartet have to be considered. This 
includes the production of resonances like K^*, N* and A++(1232). A sim- 
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plified model, assuming only an incoherent sum of these contributions is 
finally used to describe the experimental data and extract the acceptance 
corrections. This model reproduces the experimental data for many kine- 
matical variables. 

The analyzed data were taken with the High Acceptance Di- Electron Spectrometer 
(HADES) 9|] at GSI in Darmstadt, Germany. In this beam time a proton 
beam of 3.5 GeV kinetic energy was incident on a liquid hydrogen target and 
a total statistic of about 1.2 billion events was collected. 

2. Data analysis 

2.1. Evaluation of resonances contributions 

The presented analysis concentrates on the production of the A(1405) 
together with a proton and a followed by the decay of the A (1405) into 

p + p^ A(1405) +P + K+ ^ (S^vr^) + p + ^ ((n7r±)7r^) + p + (1) 
The general analysis steps to extract the A(1405) signal are presented in 



detail in 



IG 



ll| . These steps consist in identifying the four charged final 
state particles {p, K^, 7r+, tt^) and the reconstruction of the neutron via the 
missing mass to the four particles. The neutron component can be enhanced 
by an appropriate cut on the corresponding missing mass. After this se- 
lection, the S"*" and S~ hyperons are reconstructed via the missing mass of 
p, K~^, IT' or p, K~^, vr"*", respectively. By extracting the hyperon signals in the 
two spectra, the data sample is further purified, and at the same time it is 
divided into two subsamples. One subsample consists mainly of events with 
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an intermediate S"*" hyperon, whereas the other subsample contains mainly 
events with an intermediate S~ signal. For these two samples the missing 
mass spectrum of the proton and the {MM{p, K~^)), where the A(1405) 
is expected to show up, is studied separately, see fig. [TJ To show the pure 
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Figure 1: (Color online) Missing mass distribution of proton and for the two different 
subsamples within the HADES acceptance. Panel a) for events showing an intermediate 
and panel b) for events showing an intermediate E^. Experimental data (black dots) 
are compared to simulations. See text for details. 



physical signal, in both pictures the misidentification background is already 
subtracted. The treatment of this misidentification background is discussed 
extensively in 12|. The data (black dots) are compared to a sum of simu- 
lations. The strengths of the different contributions were determined by a 
simultaneous fit to four different observables, namely the two missing masses 
in fig. [Hand the two p, K~^, vr^ missing mass spectra where the missing mass 
of and S~ are visible. Details about the fitting procedure to the four 
spectra can be found in To give a full description of the experimental 
data, several contributions have to be taken into account in the simulation. A 
list of considered channels with particles in the same final state as in reaction 
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([T]) (p, ^ vr"*", vr ) is given in tabled] The channels 9 and 10 can be rejected 



Channel 


p + p ^ 


uategory 


1 


A(1405) + P + (S=^7r^) +p + 




2 


S(1385)" +p + (S=^7r+) + p + 


Ltt resonant 


O 






4 






5 


S+ + fr+ + A°(1232)/iV(1440) ^ S+ + fr+ + (pvr-) 


S+TT^ non-resonant 


6 






7 


+ 7r+ + p + 




8 


S- + + A++(1232) ^ S- + + (p7r+) 


S^7r+ non-resonant 


9 


A + 7r+ + n + 




10 


-ft's + p + n + fr+ 





Table 2: Reactions taken into account for the analysis. The channels are classified into 
two main categories. See text for details. 

n 

from the data sample, as demonstrated in jlOj]. The other channels, however, 
contain all the same final and intermediate state particles and can therefore 
contribute to the data in fig. [H They are classified into two categories: 

• "Ett resonant" are all channels, where the S and the tt stem from the 
same mother particle. They should be visible as resonances in both 
MM(p, spectra of fig. [H 

• "S+vr" (S~7r"'") non resonant" are channels which have a S+vr" (S"??"*") 
pair as an intermediate state, but the two particles are not stemming 



6 



from a common mother particle. These channels give a broad, phase 
space like distribution in the spectra of fig. [1] The "S^'^vr^ non reso- 
nant" channels can only contribute to fig. [T]a), whereas the "S^7r+ non 
resonant" channels can only give significant contribution to fig. [Ub). 

Indeed, clear peak structures around 1400 MeV/c^ and 1500 MeV/c^ can 
be observed in both spectra of fig. [H They are attributed to the channels 
1-3. The A(1520) (green histograms) is well separated from the A(1405) 
mass area and can therefore be isolated. However, the 2(1385)*^ (violet his- 
tograms) overlaps completely with this area, and it is impossible to separate 
the A(1405) and the S(1385)° in this data sample. Only in the neutral decay 
channels the two resonances show different properties (A(1405) — > S'^tt'^, BR 
33.3% and S(1385)° Avr^, BR 88%). This allows to analyze them inde- 



pendently, 
reported in 



he obtained results for the p+p data at E^-in = 3.5 GeV are 
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14l ] and yield a cross section ratio of crA(i405)/crE(i385)o ~ 1- 
This value is used as an external constraint for the analysis presented here. 
It gives the contributions of A(1405) and S(1385)° shown in fig. [1] For the 
simulation of the A (1405) a Breit-Wigner distribution (black histograms) 
was used. However, for a good agreement between simulation and experi- 
ment, the Breit-Wigner had to be simulated with a pole mass of around 1385 
MeV/c^, which then results in the solid gray histograms. 
To describe the spectra in fig. [T] completely, also phase space like distribu- 
tions (red histograms), coming from the "Svr non-resonant" channels, are 
needed. A priori it is not clear to which extent the different channels in ta- 
ble [T] contribute, as the spectra in fig. [T]are not sensitive to this information. 
Therefore, other observables have to be studied. 
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Fig. concentrates on the subsample with an intermediate S"*" hyperon. It 
shows the same data set as in fig. [T]a), but before subtracting the misidenti- 
fication background (blue histograms). Fig. [2] a) displays the invariant mass 
distribution of the proton and the vr" {M{p, 7r~)), where possible A and A^* 
resonances should appear. For extracting a possible contribution of a K*^, 
the invariant mass distribution of the and the vr" {M{K^, vr~)) is shown 
in fig. [2]b). Compared to the experimental data are simulations, where the 

150| 
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data 
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- Zitnon. res. 
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Figure 2: (Color online) a) Invariant mass of proton and tt~ and b) invariant mass of 
and vr^ for the subsample with an intermediate hyperon. The spectra show the results 
within the HADES acceptance. 

scaling of the different channels is known from the simultaneous fit men- 
tioned above. For the "S+vr^ non-resonant" contribution only channel 4 is 
included. This assumption gives already a rather good description of the 
data. As indications neither of A/N* nor of K*^ resonances are visible, only 
this channel 4 is used in the further analysis. Indeed, also the "S+tt^ non- 
resonant" part shown in fig. [1] contains only this channel. However, possible 
contributions due to the channels 5 and 6 can not be excluded completely 
by this analysis. For example the production of a K*^{892) via channel 6 
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is only slightly above threshold and thus the cross section might be just too 
small to see a clear contribution to fig. [Hb). 

To identify the different contributions to the "S^7r+ non-resonant" part, the 
invariant mass of the proton and the n'^ {M{p, 7r+)) is studied in fig. [31 Only 




Figure 3: (Color online) Invariant mass of the proton and the 7r+ within the HADES 
acceptance. Panel a): only channel 7 is used for the E~7r+ non- resonant part. Panel b): 
according to a minimization, only channel 8 is used for the S^7r+ non-resonant part. 

the data subsample of fig. [T]b) with an intermediate S~ hyperon is investi- 
gated. The misidentification background is not subtracted. Panel a) of fig. [3] 
shows the result where the "non resonant" simulations contain only channel 
7. The scaling factors for the different channels are again known from the 
fitting procedure. The data show an enhanced structure, which can not be 
described by the simulations. As a comparison, fig. E] b) shows exactly the 
same data, but now including channel 7 as well as channel 8 into the simu- 
lations. The relative contribution of these two channels is a free parameter, 
which is obtained by a fit to the experimental data points in fig. 121 The 
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fit results in a negligible contribution of channel 7. With the inclusion of the 
A++(1232) the experimental data can be described rather well. Due to this 
result, it is concluded that the "Z1^7r+ non-resonant" contribution is domi- 
nated by A'^+(1232) production. Therefore, only the channel 8 is used in the 
simulation, which is already taken into account in fig. [Hb). 
With the presented analysis a simulation model with several contributions is 
obtained, which gives reliable descriptions of the observables investigated so 
far. However, the goal of this analysis is to understand and to describe the 
experimental data within the full HADES acceptance. This asks for accep- 
tance corrections. For this purpose it might be not sufficient to study only 
invariant mass distributions, as they are not very sensitive to e.g. angular 
distributions of the produced particles. Therefore, angular distributions in 
the Center-Mass System (CMS), Gottfried- Jackson system and helicity sys- 
tem are studied in the next part. Detailed information about the properties 
of these frames and the corresponding angular distributions can be found in 
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2.2. Angular distributions 

Starting point is again reaction ([T]). Here, three particles are produced in 
the entrance channel (A(1405), p and K^). The momentum of the possible 
A(1405) is reconstructed via the missing four-vector to the proton and the 
{MV{p, K)). It is clear from the results above that this hypothetical 
particle does not always refer to a A (1405), but can also stem from all other 
channels of table [TJ Fig. H] and [5] show all angles between the three momenta 
in the three different frames for the subsample with an intermediate S"*" or 
S~ hyperon, respectively. The nomenclature is the following: 

10 




Figure 4: (Color online) Angular distributions within the HADES acceptance for events 
with an intermediate S+ hyperon (top raw: distribution of MV{p, K), p and in the 
CMS, middle raw: helicity angles of MV{p, K), p and K^, bottom raw: Gottfried- Jackson 
angles of MV{p, K), p and K+ .) 

• ^CMS- s-iigls between particle A and the beam(target) direction in the 
Center-Mass System. 

• O^Zc'- angle between particle A and particle B in the reference frame 
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Figure 5: (Color online) Angular distributions within the HADES acceptance for events 
with an intermediate hyperon (top raw: distribution of MV{p, K), p and in the 
CMS, middle raw: helicity angles of MV{p, K), p and K^, bottom raw: Gottfried- Jackson 
angles of MV{p, K), p and K+.) 

where particle A and particle C are going back to back and have equal 
momenta (helicity angle). 

• ^A-c- s-iigls between particle A and the beam/target-type protons (bt) 
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in the reference frame where particle A and particle C are going back 
to back and have equal momenta (Gottfried- Jackson angle). 

By permutation of particles A,B and C nine different observables are ob- 
tained, where some of them are not kinematically independent. The dif- 
ferent panels in fig. H] and O show the comparison between experimental 
data and simulations. A reasonable agreement could not be obtained by us- 
ing only phase space simulations for the different channels of tab. [1] The 
production of the S(1385)"^ in the CMS was found to be rather anisotropic 



12|. The production of the E(1385)° is assumed to show the same behavior. 
Therefore, the simulation of channel 2 was folded with an angular distri- 
bution in cos (^^cMs''^^^ ■ Additionally, the data sets in fig. HJa) and Oa) 
were subdivided into several angular regions of 0^^^^g'^\ These subsamples 
were analyzed independently. In this way it was found that the A (1405) and 
A(1520) seem to be produced rather isotropically in cos [OcM^t^^^ ^ whereas 
the "S+7r^ (S^7r+) non resonant" channels, namely channel 4 and 8, show an 
anisotropic behavior. This anisotropy is included by folding the simulation of 
these channels in cos {Ocm^s '^^^ with the corresponding distributions. The 
resulting distributions are included in fig. H] and |5l A reasonable agreement 
in all 18 pictures is obtained. 



3. Summary 

We studied the production of the Tr^-K^pK^ particle quartets in p+p 
reactions. It was possible to describe the experimental missing mass distri- 
bution of the proton and by a sum of different simulations, including 
A(1405), S(1385)°, A(1520) and "non-resonant Stt" production. Studying 
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several invariant mass spectra, the "non-resonant E+7r~" production seems 
to stem mainly from the reaction p + p E"^ -|-7r~ + p + K'^ and no clear in- 
dication of intermediate resonances Uke A^*/A° (p + p — ^ A^*/A° + S+ + 
or K^* {p + p ^ E+ + p + K*'^) could be seen. However, the "non-resonant 
E^TT"*"" production turned out to come exclusively from the reaction — >■ 
A++(1232) -|-E~ Furthermore, by including an anisotropic production 

of the E(1385)° channel as well as of the "non-resonant Ett" channels, our 
simulations can describe the measured data for several angular distributions. 
The overall good agreement between experimental data and simulations is a 
necessary precondition to use the obtained simulation model for acceptance 
and efficiency corrections and to finally extract cross sections for the different 
channels. 
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